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Abstract 

Transformation theory is developed for the equations of hnear anisotropic 
elasticity. The transformed equations correspond to non-unique material prop- 
erties that can be varied for a given transformation by selection of the matrix 
relating displacements in the two descriptions. This gauge matrix can be chosen 
to make the transformed density isotropic for any transformation although the 
stress in the transformed material is not generally symmetric. Symmetric stress 
is obtained only if the gauge matrix is identical to the transformation matrix, in 
agreement with Milton et al. [1]. The elastic transformation theory is applied 
to the case of cylindrical anisotropy. The equations of motion for the trans- 
formed material with isotropic density are expressed in Stroh format, suitable 
for modeling cylindrical elastic cloaking. It is shown that there is a preferred 
approximate material with symmetric stress that could be a useful candidate for 
making cylindrical elastic cloaking devices. 
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1 Introduction 



Interest in cloaking of wave motion has surged with the demonstration of the possibihty 
of practical electromagnetic wave cloaking [2] . The principle underlying the effect is the 
so-called transformation or change-of- variables method [3l H] in which the material pa- 
rameters in the physical domain are defined by a spatial transformation. The concept 
of material properties defined by transformation is not restricted to electromagnetism, 
and has stimulated interest in applying the same method to other wave fields. The first 
applications in acoustics were obtained by direct analogy with the electromagnetic case 
[3 El [7] . It was quickly realized that the fundamental mathematical identity behind 
the acoustic transformation is the equivalence |8] of the Laplacian in the original coor- 
dinates to a differential operator in the transformed (physical) coordinates, according 
to DivGrad/ — > Jdiv J~^Fi^*grad/, where F = dx/dX is the deformation gradient 
of the transformation (see Q and J = deti^. This connection, plus the realization 
that the tensor within the operator can be interpreted as a tensor of inertia means 
that the homogeneous acoustic wave equation can be mapped to the equation for an 
inhomogeneous fluid with anisotropic density. 

However, the material properties for acoustic cloaking do not have to be identified 
as a fluid with a single bulk modulus and a tensorial inertia. There is a large degree 
of freedom in the choice of the cloaking material properties [HI [lO]: a compressible 
fluid with anisotropic density is a special case of pentamode materials [TTl [T2] with 
anisotropic inertia. The non-uniqueness of the material properties (for a given trans- 
formation) is a feature not found in the electromagnetic case, where, for instance, the 
tensors of electric permittivity and magnetic permeability are necessarily proportion- 
ate for a transformation of the vacuum. The extra freedom in the acoustic case means 
that either or both of the scalar parameters, density and elastic stiffness (bulk modu- 
lus), can become tensorial quantities after the transformation. While most papers on 
acoustic cloaking have considered materials with scalar stiffness and tensorial inertia, 
e.g. O El [71 [131 [3 [13 [HI [HI [18], see [19] for a review, there is no physical reason for 
such restricted material properties. Cloaking with such materials requires very large 
total mass [9l [20] , but the more general class of pentamode materials with anisotropic 
inertia does not have this constraint. In fact, it is often possible to choose the material 
properties so that the inertia is isotropic, in which case the total mass is simply the 
mass of the equivalent undeformed region [10]. This property can be realized if the 
transformation is a pure stretch, as is the case when there is radial symmetry. This 
distinction between the cloaking material properties is critical but, judging by the con- 
tinued emphasis on anisotropic inertia in the literature, does not seem to have been 
fully appreciated. Apart from [HI [10] there have been few studies [211 [22] of acoustic 
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cloaking with anisotropic stiffness. 

The non-uniqueness of the transformed material properties found in the acous- 
tic theory transfers to elastodynamics. The first study of the transformation of the 
elastodynamic equations by Milton et al. [1] concluded that the appropriate class of 
constitutive equations for the transformed material are the Willis equations for material 
response. The general form of the Willis equations are [H |23] 

dwa=p, (la) 
a = C'^*e + S'^* u, (lb) 
p = S'^^ *e + p'^*u, (Ic) 

where e = |(Vm + (Vm)*), * denotes time convolution and is the adjoint. The stress in 
( llbp is symmetric, and the elastic moduli enjoy all of the symmetries for normal elastic- 
ity, viz. C^f^j = C^f^i and C'^f^i = C°j|;. Brun et al. [21] considered the transformation 
of isotropic elasticity in cylindrical coordinates for the particular transformation func- 
tion used by [31 H] and found that the transformed material properties are those of a 
material with isotropic inertia and elastic behavior of Cosserat type. The governing 
equations for Cosserat elastic materials [25] are 

divo- = p^^u, (2a) 
a = C'^'Vu, (2b) 

with elastic moduli satisfying the major symmetry C^f^j = C^^/^i but not the minor 
symmetry, Cjf^i ^ Cij^ki- This implies that the stress is not necessarily symmetric, 
a* 7^ a, and that it depends not only on the strain e (the symmetric part of the 
displacement gradient) but also upon the local rotation u = |(Vw — (Vu)*). Not only 
are the parameters such as C^^ in eqs. (JT]) and ([2]) different, the constitutive theories are 
mutually incompatible: one has symmetric stress, the other a non-symmetric stress. 
We show in this paper that both theories are possible versions of the transformed 
elastodynamic equations, and that they are only two from a spectrum of possible 
constitutive theories. Apart from the two references mentioned [H |2l] , the only other 
example of transformation elasticity concerns fiexural waves obeying the biharmonic 
equation [26], which is beyond the realm of the present paper. 

The purpose here is to consider the transformation method for elastodynamics, and 
to describe the range of constitutive theories possible. The starting point is the obser- 
vation [To] that the extra degrees of freedom noted for the acoustic transformation can 
be ascribed to the linear relation between the displacement fields in the two coordinate 
systems. This "gauge" transformation introduces a second matrix or tensor, in addi- 
tion to the deformation gradient from the change of coordinates. As discussed in [10] , 
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the variety of acoustically transformed material properties arises from the freedom in 
the displacement gauge. The same freedom is also present in the elastic case, and as 
we will show, it allows one to derive a broader class of constitutive properties than 
those suggested by Milton et al. [1] and by Brun et al. [24]. The material properties 
found in these studies correspond to specific choices of the gauge matrix. 

Cloaking is achieved with transformations that deform a region in such a way that 
the mapping is one-to-one everywhere except at a single point, which is mapped into the 
cloak inner boundary. This is a singular transformation, and in practice, the mapped 
region would be of finite size, e.g. a small sphere, for which the mapping is everywhere 
regular. Our objective here is to understand the nature of the material necessary to 
produce the transformation effect, in particular, what type of constitutive behavior is 
necessary: such as isotropic or anisotropic inertia. 

The outline of the paper is as follows. The notation and setup of the problem 
are given in ^ where the displacement gauge matrix is introduced. The general form 
of the transformed equations of motion are presented in ^ Constitutive equations 
resulting from specific forms of the gauge matrix are discussed in §U particularly the 
Willis equations and Cosserat elasticity, which are shown to coincide under certain 
circumstances. The special case of transformed acoustic materials is discussed in ^ 
The elastic transformation theory with isotropic density is applied in ^ to radial 
transformation of cylindrically anisotropic solids. Based on this formulation, an elastic 
material with isotropic density and standard stress-strain relations is proposed in §6.31 
as an approximation to the transformed material. A summary of the main results is 
given in ||71 

2 Notation and setup 

Two related configurations are considered: the original Q, and the transformed region 
CO, also called the physical or current domain. The transformation from f2 to w is 
described by the point-wise deformation from X & Q to x & u. The symbols V, Vx 
and div, Div indicate the gradient and divergence operators in x and X, respectively, 
and the superscript t denotes transpose. The component form of divE is dEi/dxi 
or dEij/dxi when E is a vector and a second order tensor-like quantity, respectively. 
Upper and lower case subscripts {I , J, . . . ,i, j, . . .) are used to distinguish between 
the domains, and the summation convention on repeated subscripts is assumed. It is 
useful to describe the transformation using language of finite deformation in continuum 
mechanics. Thus, X describes a particle position in the Lagrangian or undeformed 
configuration, and x is particle location in the Eulerian or deformed physical state. 
The transformation or mapping is assumed to be one-to-one and invertible. For perfect 
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cloaking the transformation is one-to-many at the single point X = but this can 
be avoided by always considering near-cloaks, where, for instance, the single point is 
replaced by a small hole which is mapped to a much larger hole. 

The deformation gradient is defined F = V xX with inverse F~'^ = VX, or in 
component form Fu = dxi/dXj, Ff-^ = dXj/dxi. The Jacobian of the deformation 
is J = deti^ = \F\, or in terms of volume elements in the two configurations, J = 
dv/ dV. The polar decomposition implies F = VR, where R is proper orthogonal 
= = /, deti? = 1) and the left stretch tensor V E Sym"^ is the positive 
definite solution of = B where B (= FF^) is the left Cauchy-Green or Finger 
deformation tensor. 

The infinitesimal displacement U{X,t) and stress E(X,t) satisfy the equations of 
linear elasticity in the original domain: 



DivS = pqU, 



in n, (3) 



where po is the (scalar) mass density and the the elements of the elastic stiffness tensor 
satisfy the full symmetries 

^(0) _ ^(0) ^(0) _ ^(0) / 

^IJKL — ^JIKL^ ^IJKL — ^KLIJ- V^) 

The first identity expresses the symmetry of the stress and the second is the consequence 
of an assumed strain energy density function. The total energy density is the sum of 
the strain and kinetic energy densities, 

£o = m + % ^\iemWo = \Cf]j,^UjjUL,K, % = \poU'u. (5) 

Particle displacement in the transformed domain is u{x, t). Our objective is to find 
its governing equations. In order to proceed, we need in addition to the geometrical 
quantity F , a kinematic relation that relates the displacements in the two domains. We 
assume a linear "gauge" change in the displacement defined by a non-singular matrix 
A as 

U = A'u {Ui = AuUi). (6) 

According to its definition the matrix A is, like F, not a second order tensor because it 
has one "leg" in both domains. The choice of the transpose, j4* in equation (|6]), means 
that A and F are similar objects, although at this stage they are not related. 

The arbitrariness in the choice of A is the central theme of this paper. This approach 
generalizes that of [10] which was restricted to acoustic materials, and of Milton et 
al. [1] for elasticity. The point of departure with [1] here rests with the assumed 
independence of the gauge matrix A. Milton et al. assume a similar relation between 
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the displacement fields; eq. is identical to eq. (2.2) in [1]; however, the matrix 
^ in [1] is assumed at the outset to be equal to the deformation gradient (A = F). 
We will return to this distinction later. As noted by Milton et al. [1], the relation 
dX = F^^ dx might lead one to expect A = F"*" by identifying dX and dx with U 
and It, respectively. However, the displacements are not associated with the coordinate 
transformation, unlike in the theory of finite deformation, and hence F and A are 
independent. Milton et al. [T] specify A = F on the basis that this is the only choice 
that guarantees a symmetric stress. We will return to this point later. 

3 General form of the transformed equations 

Under the transformation and the gauge change the energy density transforms as £q — )• 
£ = W + T according to £^ d V = d Vq, so that 

8 = W + r= I {CJjA'L (^i^jj) ,i {-^i^il) ^k^iiFkK + poiiiUjAiAjj. } (7) 

Hence, 

W = iJ-'CfJ^L FuFuK {ujA^j) [uiAil] T = \u*fm, (8) 

where the (symmetric) density tensor is 

p = p' = p,J-'AA\ (9) 

The equations of motion in the deformed, or current material, are determined by the 
Euler-Lagrange equations of the Lagrangian density £ = W — T, as 

A, J {j-^Cf]j^L Fil{uiAiL),kFkK)^^ - P^jUi = 0. (10) 

Using the identity 

{J-'Fu), = 0, (11) 

this can be written 

QijIJ {JC^JkL {uiQklKL),k) i - PijUi = 0, (12) 

where the fourth order quantity 

QijIJ = J"^FuAjj, (13) 

is introduced for later use. 

The transformed system ( fT2l) is equivalent to the equilibrium equations 

(Jij^i=Pj, (14a) 
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and the constitutive relations 

^ij = Clfkini,k + Stfi ill, pi = Stfi u,,i + pfiik, (14b) 
with parameters defined as follows in the Fourier time domain (dependence e~^^) 

Qfkl = JCfjKLQijIjQklKL, (15a) 

St^ = i-tcuy' JCi%^Q,,jjQ,iKL,k, (15b) 

Pjf = Pjl + {-i^y^ JCf]KLQijIJ,iQklKL,k, (15c) 

where the density pji is given by ([9]). The elastic moduli and the density satisfy the 
general symmetries 

^ifki = '^mj^ Pjf = Pif' (16) 

but not the full symmetries required for the Willis constitutive model ([1]). Equations 
(IT^ - fllSp are the fundamental result of the transformation theory. The remainder of 
the paper is concerned with their simplification and interpretation. 

Note that the transformed stiffness may be expressed in a form similar to the Kelvin 
representation for the tensor of elastic moduli [27] , as 

6 

C^ff = ^i^(")s(°) ® S(°\ = JK'^^\ S(") = J-iFP(")A*, (17) 

where k'^'^ > are the Kelvin moduli, P*-"^ G Sym, tr P*-°^P^'^-' = Saf3, are the eigen- 
strains/eigenstresses, such that the original stiffness has the unique decomposition 

6 

c(°) = j];i^i")p(")®p("). (18) 

The transformed matrices S*-"^ do not inherit the orthogonality of the original eigen- 
strains/eigenstresses P'-"'*, so that f|T71) is not the exact Kelvin representation in the 
transformed coordinates. It does, however, illustrate that the transformed stiffness is 
positive definite, even though S*-"-* are in general not symmetric. The representation 
(fT7|l is particularly useful in the limiting case of an acoustic fiuid in the original domain 
for which only one of the K^"^ is non-zero, discussed later. 

4 Transformed equations in specific forms 
4.1 Willis equations: A = F 

The absence of the minor symmetries under the interchange of i and j in C^^/^i and 
Sifi of (fT5|) implies that the stress is generally not symmetric. Symmetric stress is 
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guaranteed if Qijij = Qjuj (see eq. (fT3l) ). which occurs if the gauge matrix is of the 
form A = (F, for any scalar C 7^ 0, which may be set to unity with no loss in generality. 
This recovers the results of Milton et al. [1] that the transformed material is of the 
Willis form, eq. ([T]). As noted in [1], this is the only choice for A that yields symmetric 
stress. 

In summary, the governing equations are (HM with material parameters defined by 
(Hg) and 

Q^JIJ = J-'FuFjj. (19) 
The parameters now display the full symmetries expected of a Willis material: 

/^cS /^cff /^cff /^cS eff off Qcff Qeff (Of\\ 

^ijkl — ^klij^ ^ijkl — ^jikl^ Pjl — Plj '-'ijl — '-'jil- ■ \^^) 

Note that the stiffness tensor is 

= FuFjjFkKFiL J-'cf]^L = VuVjMkVil J-^Cf]^^, where 

^fjKL = RimRjnRpkRlQ Cfjj^pQ (21) 

are the original moduli in the rotated frame. The full symmetry of the stiffness 
tensor also follows immediately from the representation f|T71) with symmetric 8^°"^ = 

4.2 Cosserat elasticity: A = constant 
4.2.1 General form 

The constitutive parameters ( ITSjl simplify considerably if the fourth order quantity Q 
satisfies 

Qi,ij,i = 0. (22) 

This differential constraint combined with (fTTl) implies that the gauge A must be con- 
stant. In that case the transformed equations of motion become 

= Pij ^jj (^ij = d^jj^i Ui^ki (23) 

where the effective elastic moduli are defined by (fT3l) and (11 Sal) and the density tensor 
p is given in ([9]). 

Note that the elastic moduli satisfy the symmetry (fT6l) i associated with the trans- 
formed energy density W = \C^f^.iUj^iUi^k- But Cff^^i does not satisfy the minor symme- 
try (j4j)i since 



1^(0) 

J'^Cf]^^ FukAil {FuA.j - FjjAu) , (24) 



Qfkl - Cjikl - J ^ClJKL^kKAiL{FiiAjj - FjjAij) 

-ir^(o) 



8 



which is non-zero in general (note that the second form in ( 12^ uses the minor symmetry 
for the original moduli C}j^^). This means that the stress is not necessarily 
symmetric, cr 7^ cr*, which places the material in the framework of Cosserat elasticity 
[25]. The number of independent elastic stiffness elements is at most 9(9 + l)/2 = 45 
as compared with 6(6 + l)/2 = 21 for normal linear elasticity. 

4.2.2 Cosserat elasticity with isotropic density: A = I 

Isotropic density can be achieved by taking the constant matrix A proportional to the 
identity, A = (I, with C = 1 without loss of generality. In this important case we have 
p = pi, with 

P = PqJ \ Qjfci = J~^Cjfj^iFiiFkK- (25) 

4.3 Examples 

4.3.1 Example 1: SH motion in a plane of material symmetry 

The original moduli are assumed to have a plane of symmetry perpendicular to the 
X3-axis, and the transformation is assumed to preserve the out of plane coordinate: 
X3 = X3. Consider shear horizontal motion U = (0, 0, t/(Xi, X2, t)) satisfying the 
scalar equation 

{C'^LUb)^^ = PoU, (26) 

with indices A,B^ {li2}. Under these circumstances, the equation for SH motion 
in the transformed domain, u = {0,0,u{xi,X2,t)), is the same for both the Willis 
constitutive equations {A = F) and the Cosserat model with isotropic density {A = I). 
Thus, 

(C:3/33«,/3),, = P^, (27) 

where a,/? G {1,2}, p = J^^po, and C^|^3 = J^^C^^l^^FaAFisB- The equivalence may 
be expected since the only relevant element of the gauge matrix, ^33, is the same for 
both models (A33 = 1). 

The above conclusion for SH motion in the presence of orthotropic moduli relies 
only upon the scalar nature of the motion in the original and transformed domains. As 
such, the SH results also follow from those in ^ for fluid acoustics under the standard 
replacement of fluid density and bulk modulus with inverse shear modulus and inverse 
solid density, respectively. 

4.3.2 Example 2: F = constant 

When both A and F are constant the Willis equations simplify to those of normal linear 
elasticity {Siji = 0). At the same time, the constant deformation gradient implies that 
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A = F (=constant) is a permissible choice for the Cosserat medium. The Wilhs and 
Cosserat materials are then coincident with density 

p = p' = p^.r^B. (28) 

and fully symmetric elastic moduli given by fl2T]) . Note that the density p is anisotropic 
unless B = al, which means the deformation is a pure expansion, possibly with rota- 
tion. But this is a rather trivial case. 

Consider an isotropic initial material with original moduli C}j^^ = Xq^ij^kl + 
fJ'oiSixSjL + SilSjk), or equivalently, 

C(°) = Ao/®/ + 2/io/K/, where {X^X)Y = ^X{Y + Y')X. (29) 

The rotated moduli of fl2Tl) 9 are therefore unchanged, C^^^ = C^'^\ and the current 
density and moduli are 

C = \B ®B + 2fiB^B, p = pB, where {A,/i,p} = J"^ {Ao,/io,po}- (30) 

It is of interest to consider plane wave motion, u = qg(n^x — vt) for unit vector n, 
constant q, and g E C^. The equation of motion ( !23l) i implies that the polarization 
vector q satisfies 

{Q - v^p) q = 0, where Q = {X + fi) {Bn) ® (Bn) + fi{n*Bn) B. (31) 
The solutions of fl3T]) are of two distinct types: 

longitudinah g || n, v = {n^BnY^^ cl, ^1 = ^^ = ^^^^, (32a) 

P Po 

quasi-transverse: q _L Bn, v = (n^Bn)^^'^ ct, ci = — = —. (32b) 

P Po 

The slowness vector is s = f ~^n, and the slowness surface is the envelope of the 
slowness vectors for all propagation directions n. The slowness surface is therefore 
comprised of two sheets which are similar ellipsoids: s^Bs = c~^, a = L,T. Note that 
the polarizations do not in general form an orthogonal triad. 

5 Acoustics as a special case 

5.1 General form of transformed equations 

We consider the simpler but special case of acoustics in order to further understand 
the structure of the elastic transformation theory. The acoustic equations are unique 
in the sense that they are the single example of a pentamode material [TTl [12] com- 
monly encountered in mechanics. We will demonstrate that the pentamode property 
introduces a unique degree of freedom not available in the fully elastic situation. 
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The elastic stiffness of an acoustic fluid is 

Co = i^o I ® I, (33) 

which is of pentamode form [TT], i.e. the 6x6 Voigt matrix associated with the elements 
^fjKL ~ Ko6ij6kl has five zero eigenvalues. Equation f|T2|) becomes, using fl33|) . 

Qij [JKq {uiQki),k) ^ - pijUi = 0, (34) 

with the density tensor given by , and 

Q = J-^FA* (g,, = J~'FuAji = QijKK) ■ (35) 

The general form of the governing equations are again of the form (fT4|) . with material 
parameters 

C^ifki = KoJ^^BikAjNAiN, (36a) 
S^fi = {-iujy'KoJ-'BikAjNAiN,k, (36b) 
Pjf = Pji + {-iu})~^KoJ^^BikAjN,iAiN,k, (36c) 

recalling that B = FF^ . The stress is again not generally symmetric, unless A is 
proportional to F, in which case the transformed material is of Willis form, see §4.11 



5.2 Cosserat and related forms of the transformed equations 

The simplified form of the elasticity in the acoustic fluid implies that the condition (|22|) . 
which is required to simplify the form of the transformed equations, itself simplifies to 
the condition that Q be divergence free: 

divQ = (g^,, = 0). (37) 

Note that this is a necessary but not sufficient condition for the more general elastic- 
ity version f l^ . Assuming that holds, the transformed equations flM]) have the 
simplified Cosserat structure ( 12^ . with 

p = Po^Q*B"^Q, = KqJCi®Ci, divQ = 0. (38) 

We note the symmetries p = , C^ff^i = Clf^j, but the minor symmetry = Cjfi^i 
does not in general hold unless Q is symmetric [TD]. Thus, the transformed acoustic 
equations are those of a pentamode material of Cosserat type with anisotropic density. 
All previous studies of transformation acoustics assumed a priori that the transformed 
materials must have symmetric stress. The present results show that the more general 
structure of the transformed properties is that of a material with stress not necessarily 
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symmetric. The pentamode structure of C° implies that the equations of motion (123|) 
can be expressed in a form that is clearly related to acoustics [9] by using a scalar 
"pseudo-pressure" p and "bulk modulus" K = KqJ, 

pw = -QVp, p = -frtr(QVv). (39) 

The condition f l37p can be achieved, as in the elastic case, with constant gauge A. 
For instance, taking A = I, yields material properties (see (l25l) ) 

p = poJ~\ C'^ = KqJ-^ F ® F. (40) 

This describes a material with isotropic density of general pentamode/Cosserat form. 
That is, the stiffness is pentamode (a single nonzero eigenstiffness |27j) and a sin- 
gle eigenstress of generally non-symmetric form, hence Cosserat. As discussed in [5], 
isotropic density with symmetric stress can be achieved if div KV = for some func- 
tion /i(x). One important case is when the deformation is a pure stretch F = V (with 
h = J~^, see eq. (fTTj) ). for which the transformed material is pure pentamode with 
isotropic density: 

p = J-i I, C = KqJ'^ V (g) V, pure stretch. (41) 

More general conditions under which pure pentamode material with isotropic density 
can be achieved are discussed in |10j . 

Condition (1371) may also be satisfied by non-constant gauge matrices. For instance, 
A = JF~* gives Q = I which clearly satisfies fl37|) . In this case the transformed medium 
has the properties 

p = PoJB-\ C'^ = KoJl®I. (42) 

This corresponds to a fluid with isotropic (hydrostatic) stress, bulk modulus K = KqJ, 
and anisotropic density p. This type of material was the first to be considered for 
acoustic cloaking [5], [6l [7] yet it is remarkable that it does not possess a generalization 
to elasticity. Rather, it is made possible by the simple structure of the second order 
quantity Q as compared with its fourth order elasticity analog Q. The condition (12^ 
for Q is only satisfied by constant A, but the analogous condition fl37|) for Q has at 
least one non-constant solution for A. The one noted here, A = JF~*, yields the class 
of acoustic cloaking materials that use anisotropic inertia as the active mechanism. 

The transformed acoustic material can also be understood as the special case in 
which five of the Kelvin moduli K^"'^ vanish (see (fTSjl ). say a = 2,3, ... ,6, and the 
remaining single eigenvector is the identity, P*-^^ = I. The transformed material has a 
single non-zero eigen-stiffnesses with associated eigen-tensor (see eq. (135|) ) 

S(i) = j-iFA*. (43) 

The inertial transformed material (jUJ results from the choice of A that makes S^^-* = I, 
that is A = JF"*. 
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6 Applications in cylindrical elasticity 



The non-uniqueness in the form of the transformed equations of elasticity provides 
the designer of elastic cloaking devices with a wide variety of possible materials from 
which to choose. These range from materials with Willis constitutive behavior ([1]), to 
Cosserat materials (l25ll . each special cases of the general constitutive relations ( Il4l) - 
f|T5|) . The same non-uniqueness exist for acoustic cloaking, where designs based on 
inertial cloaking on the one hand, and pentamode materials on the other, reflect the 
choice of completely different material properties. One aspect common to all of these 
materials is their exotic nature. In the absence of available materials with exactly the 
right properties it is reasonable to ask what can be achieved using "normal" materi- 
als, i.e. those with isotropic density and standard linear elastic response, including a 
symmetric stress. 

With that goal in mind we examine in this section cloaking in the context of the 
Cosserat materials described by fl25|) . This class of metamaterials is chosen as the 
starting point because of its property of isotropic density and the fact that the con- 
stitutive behavior is local. That is, the stress depends on the displacement gradient 
alone, and not on displacement as in f|T^ - f|T5|) . The goal is to find if there is a material 
with symmetric stress that provides an optimal, in some sense, approximation. We 
consider transformations in cylindrical coordinates, starting with the formulation of 
the Cosserat constitutive equations in cylindrical coordinates. Thus, for the remainder 
of the paper A is assumed to be a constant, and we use C instead of C°^. 

6.1 General theory for Cosserat materials 
6.1.1 Cosserat notation in cylindrical coordinates 

The cylindrical coordinates are referred to by the indices 1, 2, 3 for r, 9, respectively. 
The usual Voigt notation, which means {1,2,3,4,5,6} = {11,22,33,23,31,12}, is 
augmented with three additional indices to describe Cosserat elasticity: {4, 5, 6} = 
{32, 13,21}, so that the elastic stiffness tensor becomes in the 9-index Voigt-Cosserat 
notation, 
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Following [28], the traction vectors tj = ti(x, t), i = r,6,z, are defined by the 
orthonormal basis vectors {e^, e^, e^} of the cylindrical coordinates {r, 6, z} according 
to tj = ejtr [i = r,6,z), where <t(x, t) is the stress, and a comma denotes partial 
differentiation. With the same basis vectors, and assuming the summation convention 
on repeated indices, the elements of stress are 0"^ = CijkiSki where e = |(Vu + Vu*) 
is the strain, Cijki = Cijkii^x) are elements of the fourth order (anisotropic) elastic 
stiffness tensor. The traction vectors are |29] 



/t,\ 
te 



/Q R P\ 

R* T S 
VP* S* Mj 



( 



u 



i(u,,+Ku) 



v 



u 



/Q R P\ 

R* T S 
VP* S* M/ 



({erer) (cr-ee) (er-ez)\ 
V (e.e,)/ 



where K = ® e^. — e,. (8> eg, and in notation similar to that of [30], the matrix (ab) 
has components {ab)^^ = aiCijkibk for arbitrary vectors a and b. The explicit form of 
the various matrices is apparent with the use of Voigt's notation 



Q 
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C15 
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C56 


C55 
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C25 
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C45 
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, C45 


C44 


C34 




VC35 


C34 


C33/ 






C12 


cu\ 




^66 


C26 


C46 1 






C25 


C45/ 



(45) 



Note that the six elements in each of the symmetric matrices Q, T, M and the nine in 
each of S, P, R, are independent. That is, there is one-to-one correspondence between 
the elements of the six matrices and the 45 independent components of the Cosserat 
elasticity tensor 



6.1.2 Cylindrically anisotropic materials 

We consider materials with no axial dependence whose density and the elasticity tensor 
may depend upon r, i.e. p = p{r) and Cijki = Cijkiir). We seek solutions in the 
form of time-harmonic cylindrical waves of azimuthal order n = 0,1,2,... and axial 
wavenumber kz, such that the displacement-traction 6- vector is of the form 



u(r, 6, z, t) 
.irtr{r, 9, z, t) 



•qyr) e 



i(n6+kzZ—Lji) 



(46) 



where 77 depends only on the radial coordinate r. Accordingly, the governing equations 
of motion reduce to an ordinary differential equation for this 6x1 vector [29] : 



^--G(r)j7(r) 
dr r 



0, 



(47) 
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with 6x6 system matrix G, where 



iG{r) 



(T - R+Q-iR - pwV^) R+Q-i, 



-Q-^P 

I [P*Q-iR - S - (P*Q-iR - S)^ + zA;,r (P*Q-ip - M)] -P'Q \ 

R = Rk, S = kS, T = k+Tk = T+, K = K + mI = -«+. 



Note the symmetry G = JG"'"J for real-valued material constants and u, kz, where 
J has block structure with zero submatrices on the diagonal and off-diagonal identity 
matrices. This hermiticity-like property has important physical consequences such as 
conservation of energy [28]. 



6.1.3 Transformation in cylindrical coordinates 



Let R = (Xf + Xl)-*^/^, r = (xl + x'^Y^'^, and consider the reverse deformation R 
R{r),X3 = X3. The deformation gradient is then 



F = alr + pie + with a = " ^' ^^^^ 



where /,. = er ® Br, le = &e ® 69, Iz = S-z® ©z- Taking the (assumed constant) gauge 
matrix as A = J gives isotropic density, p = pi and Cosserat elastic stiffness C, where 
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according to eq. ([25 
1 



P 



a/3 



Po, 



(50a) 
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12 



ac 
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13 
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14 



ac 



(0) 
15 



2^(0) 
15 
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.(0) 
-^34 
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-35 



.(0) 
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ac 



(0) 



36 



(0) 



-,(0) 
-46 



C44 



.(0) 
--45 

M 

"55 



ac 



ac 



(0) 

45 

(0) 
55 



ac 



(0) 
46 



ac 



(0) 
56 



2M oM 



a c 



55 



a c, 



S 



Y 



M 



56 



2JO) 



a/3c 



(0) 
56 



a c. 



66 



v 

Equivalently, the six matrices of 

s = -s(°), 

a 



a/3c66 



(50b) 

that are in one-to-one correspondence with C are 



^rj.(0) 



M 



a 



ap 



/3 



(51) 



These results are consistent with the system equation P7|) in the current coordinates 
and the analogous equation in the original coordinates, 



dry 



(52) 



dR R 

with G^°\R) defined in the same way as (jH]) using Q(°\ . . ., K^^\ Note that the 
transformed version of ( l52l) may be obtained directly by multiplication with the factor 
di?/dr, 



dr) 

dr 



i dR 
R~d7 



G(°)(^)^- 



Comparison of eqs. (147|) and (!53|) implies that the transformed system matrix is 

r di? 



G(r) 



dr 



G^'\R). 



(53) 



(54) 



This relation, combined with the block structure of G in ( l48i) and the analogous form 
for G^''^ implies the matrix relations ( ISTl) . 
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6.2 Example: a cylindrically orthotropic material 
6.2.1 Transformed elastic moduli 



The initial material is assumed to be cylindrically orthotropic, with stiffness in the 
usual Voigt notation ({1, 2, 3, 4, 5, 6} ^ {11, 22, 33, 23, 31, 12}) given by 



C(o) 



(0) 

SO) (0) jo) 



11 



'12 



C 



V 



13 







•-12 



-22 



23 








c 



--23 



(0) ^(0) ^(0) 



C33 











JO) 

-44 










.(0) 

-55 




o\ 






-66 / 



(55) 



Hence, using fl50bl) . the transformed elastic stiffness tensor becomes in the 9- index 
Cosserat notation, 



C= — 

a(3 



-11 



'22 



(0) 
'-33 



Y 









JO) 

'44 



M 







1-44 









SO) 

'55 
















































(0) 
«45 


















" '-66 


a/3c 



\ 



(0) 
66 
(0) 
66 ' 



(56) 



6.2.2 SH motion 

The equation of motion in the undeformed coordinates for U = (0, 0, U{r, 9, t)), 

+ ^A^uUe)^, - PoU = 0, (57) 
transforms to the following equation for u = (0, 0, u{r, 9, t)), 

^ {rciiu^r),T + i(c44^i,e),e - pii = 0, (58) 



where 



RdR _ (o)rdi? _ (o)i?di? , , 

P — Po — —, — 5 C44 — C44 — - — , C55 — C55 — - — . 
r dr it dr r dr 



This is an example of the general result in §4.3. II for SH motion in a plane of symmetry 
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6.2.3 In-plane motion 

The displacement is assumed to have the formit = (Mr(r, 6, t),U0{r, 9, t),0). The density 
is again isotropic given by eq. f l59|) i and the relevant moduli are, from eq. fl50b|) . 



Cll 



Cm 



.(0) 

-11 

.(0) 
-66 



RdR 

r dr ' 
R dr 

VJr' 



C22 



C66 



^''Rdr' 
^(0) r^cLR 

"66 D J ' 

R dr 



Cl2 



C66 



.(0) 
-12 5 

M 

-66 • 



(60a) 
(60b) 



For isotropic initial stiffness tensor, these expressions agree with the Cosserat elastic 
moduli found by [21], where only the particular transformation of [3lll] was considered, 
i.e., R = ri{r — tq) / (ri — tq) for r G (ro, ri]. 



6.3 A symmetric approximation for A;^ = 

The example in §6.21 achieves cloaking of in-plane elastic wave motion using elastic 
moduli of Cosserat form, ( l60l) . generalizing the findings of [21]. The existence of three 
distinct in-plane shear moduli, C1212, C2121, C1221 = C2112, is a property of the Cosserat 
model, but one that is difficult if not impossible to realize in practice. We now show 
that there is a preferred approximation with only a single shear modulus. 

Consider the transformed system matrix G with = for the initial cylindrically 
orthotropic material ( l55i) . which follows from ( 15^ as 



iG 



a 



/o 





4? 




/joy 



—I 



iuj'^R^PqI 



-K 









\ 


^-66 










jo)- 

•-55 


-1 

/ 





1 0^ 




•-12 





1 





0^ 





V yo 0^47 \ 0/ y 

(61) 

Let Cijki be a set of moduli in the transformed coordinates, corresponding to a normal 
elastic solid (symmetric stress) of cylindrically orthotropic symmetry, with moduli (in 
the standard Voigt notation). 



/ Cll C12 ci3 \ 

C12 C22 C23 

Ci3 C23 C33 

C44 

C55 

\ cee/ 



(62) 



Assuming the density is isotropic and equal to the transformed density of eq. fl50ap . 
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the system matrix associated with the stiffness Cijki becomes 
/ 



iG 






fen-' 








K —i 









^ 




K 



















1 0\ 


K — iuj^r'^pl — 


H 


Cll 
















0/ 



(63) 



Comparison of and ( 1^ suggests the identification 



Cll — -7.(^11 , C22 — — C22 , C12 — Cy2, 1 (^D4aj 

p a 

C44 — — C44 , C55 — — C55 , C66 — -jCgg . (^b4Dj 

a Q! p 

This set of moduh has the property that they correspond to a normal elastic material, 
as compared with the Cosserat material required for the exact solution. Comparing 
the latter, given by f lSB]) . with the proposed moduli flM|) shows that 

Cll = Cll, C22 = C22, C12 = C12, C44 = C44, C55 = C55, Cqq = Cqq. (65) 

Let ^{r) be the solution for the approximate but symmetric moduli (|62|) . 



^ - - G(r)fj(r) = 0. (66) 
dr r 

It is evident that the exact and approximate systems (l47l) and (!66|) have identical 
SH solutions. The approximate solution defined by the moduli C^ki possesses a further 
interesting property related to in-plane wave motion. Thus, the difference between the 
system matrices of the exact transformed medium and that of the approximate material 
follows as 

^G-^G = /(r)D^A, (67) 

where /W = - ^> D = (^^ _D+) ' ^i = | -m 1 0|. (6^ 

The difference is independent of the material properties, a function of only the trans- 
formation function through /(r) and of the azimuthal index n through D which is rank 
two and has the property D'^ = D or for any odd or even k. This suggests that 
the choice (!64l) of the approximate material parameters is in some sense preferred over 
others. Some consequences are discussed next. 

The matricant solutions M(r, ri) and M(r, ri) of, respectively, the exact and ap- 
proximate systems ( l47l) and ( |66l) . are defined [31] such that 77 = Mc and r] = Mc 
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where c is the arbitrary initial data vector at r = ri. The matricants therefore satisfy 
differential equations similar to the state vectors with initial conditions M(ri,ri) = I, 
M(ri, ri) = I. According to [3ll §3 of Ch.VII], the matricant ( 147|) and ( l66l) on an inter- 
val [ri,r2] are related as M = MP through the matrix P satisfying P' = M ^AMP. 
Noting that R' = AR with initial condition R(ri, ri) = I has a simple explicit solution 

R"=I + diag([(gi)*'-l]D„ [(r|l)-'_i]D*), (69) 

it follows that 

M = MRT where T' = AiT, T{n,n) = l, (70) 

and Ai = (MR)~-'^(AM — MA)R. The latter vanishes at r = ri, which suggests that 
the convergent Peano series T = Tq + Ti + T2 + . . . with Tq = I, T'j = AiT^.i, 
j = 1,2,..., forms a natural and regular perturbation solution. 

Note that in cylindrical cloaks the original domain [Ri < ri, i?2 = T2] is mapped to 
the smaller one [ri, r2] so that the elements of R in ( l69l) are sign definite. 

7 Discussion 

A complete transformation theory has been developed for elasticity. The material 
properties after transformation of the elastodynamic equations are given by eqs. (fT3|) 
through f fTSj) . The constitutive parameters depend on both the transformation and 
gauge matrices, F and A, and do not necessarily have symmetric stress. It was shown 
in ^ that a priori requiring stress to be symmetric implies that the material must 
be of Willis form ([1]), with yl = as Milton et al. [1] found. The emphasis here 
has been on exploring the consequences of relaxing the constraint of symmetric stress. 
There are several reasons for doing so. First is the fact that the transformation of the 
acoustic equation in its simplest form, i.e. by identifying an inertial tensor p = J~^FF^ 
from the differential identity DivGrad/ J div J^^FF* giad f , does not follow from 
the transformed Willis material, even though the inertial fluid has symmetric stress. 
Other types of transformed acoustic fluids are possible (e.g. pentamode materials), 
again with symmetric stress and not contained within the framework of eqs. ([1]). A 
second and more practical reason for considering the general material as defined by eqs. 
(fT3|l through (fT5|) with F and A distinct is to broaden the class of materials available 
for design of elastodynamic cloaks. 

Allowing A to be independent of leads to constitutive models that differ markedly 
from the Willis material model. In this paper we have emphasized solutions correspond- 
ing to time- independent material parameters obtained when A is assumed constant (see 
eqs. dl]), (USD and f lTSiD ). 

p = Po r^AA\ C,,ui = J-^Cf]^^ FuA.jFukAil. (71) 
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These transformed quantities correspond to a material with anisotropic density tensor 
and stress-strain relation of Cosserat type (non- symmetric stress). Setting A = I 
ensures that the transformed density is isotropic (see eq. ( l25l) ) 

P = PO, Qjkl = J'^CfjKl^ilPkK- (72) 

However, there is no general choice of the gauge matrix A that will make the stress sym- 
metric for a given transformation. This feature distinguishes the elastic transformation 
problem from the acoustic case, for which it is always possible to achieve symmetric, 
even isotropic (hydrostatic), stress. If the transformation is homogeneous, correspond- 
ing to constant F, it is possible to make the elastic stress symmetric, although at the 
price of anisotropic inertia (see eq. (l28l) ). 

Materials displaying non-symmetric stress of the type necessary to achieve elasto- 
dynamic cloaking while difficult to envisage, are not ruled out. Effective moduli with 
the major symmetry Cijki = Ckuj that do not display the minor symmetry Cijki = Cjiki 
are found in the theory of incremental motion superimposed on finite deformation [32] . 
The similarity with the transformation problem is intriguing: small-on-large motion in 
the presence of finite prestress corresponds to a transformation of an actual material via 
a deformation. The deformation in the small-on-large theory is however quite distinct 
from F in the present context. The formal equivalence of the constitutive parameters 
( 172|) with the density and moduli for incremental motion after finite prestress offers the 
possibility for achieving Cosserat elasticity of the desired form. The crucial quantity 
is the finite (hyperelastic) strain energy function of the material, which after prestress 
should have the desired Cosserat incremental moduli. Future work will examine this 
connection and the types of strain energy functions required. 

Another approach is to seek materials with normal elastic behavior that approx- 
imate, in some sense, the Cosserat material. Preliminary work in this direction has 
been considered here. The general theory has been applied to the case of cylindrical 
anisotropy for arbitrary radial transformation i? — )■ r. The equations of motion for the 
transformed Cosserat material have been expressed in Stroh format, eq. (H7|) . suitable 
for numerical implementation. The material required for cloaking of in-plane elastic 
waves is of Cosserat type with isotropic density. A normal elastic material with density 
( l5Uap and elastic moduli defined by eqs. f lB^ and flM]) appears to provide a natural 
approximation. The properties of this type of approximate material is the subject for 
planned further study, analytical and numerical. 
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